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Abstract 

We study how the entanglement of an entangled pair of particles is affected when one or both 
of the pair is uniformly accelerated, while the detector remains in an inertial frame. We find that 
the entanglement is unchanged if all degrees of freedom are considered. However, particle pairs 
are produced when a relativistic particle is accelerated, and more bipartite systems emerge, the 
entanglements of some of which may change as the acceleration. In particular, the entanglement of 
a pair of accelerating fermions is transferred preferentially to the produced antiparticles when the 
acceleration is large, and the entanglement transfer is complete when the acceleration approaches 
infinity. However, no such entanglement transfer to the antiparticles is observed for scalar particles. 
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I. INTRODUCTION 



Entanglement is an important property of quantum mechanical systems. It is useful in the 

field of quantum information and quantum computing, such as in quantum teleportation if . 

ft fl 
It also finds many applications in quantum control [2| and quantum simulations [3]. Study- 
ing quantum entanglement in relativistic systems may give us insights on the relationship 
between quantum mechanics and general relativity. It has been shown that entanglement 
is Lorentz invariant j^, 0]. However, an accelerating observer measures less entanglement 
than an inertial observer in both the scalar and fermion [8| cases. This degradation in 
entanglement is due to the splitting of space-time, as a result of which the vacuum observed 
in one frame can become excited in another frame - the case of Unruh effect [9|. Classi- 
cally, the trajectory of a uniformly accelerating particle observed by an inertial observer is 
the same as that of an inertial particle measured by a uniformly accelerating observer with 
appropriate acceleration. We are interested in how the acceleration of particles affects the 
entanglement of the originally entangled states. 

In this paper, we analyze the entanglement of accelerating particles in three cases: a 
non- relativistic wave packet, scalar and fermion particles. We compare the entanglement of 
accelerating particles, as seen 



by an accelerating detector 



jy an inertial detector, with that of inertial particles observed 
{ [h]] . We find that when all degrees of freedom are consid- 



ered, the entanglement is unchanged. However, pair production occurs when a relativistic 
particle accelerates, and there are new bipartite systems. We find that the entanglements 
of some of the new bipartite systems can depend on the acceleration. In particular, a pair 
of accelerating fermions transfer their entanglement preferentially to the produced antipar- 
ticles when the acceleration is large, and the entanglement transfer is complete when the 
acceleration approaches infinity. However, no such entanglement transfer to the antiparticles 
is observed for scalar particles. 

This paper is organized as follows. In Sec. HIl we add a potential term in the Schrodinger 
equation that would lead to uniform acceleration in the classical limit. Then we construct 
a wave packet solution and a two-body entangled wave function, and we calculate the en- 
tanglement by Schmidt decomposition as a function of acceleration. In Sec. 1111} we add 
the same potential term in the Klein-Gordon equation and Dirac equation. A wave packet 
solution can be obtained [3] to give an intuitive picture of how a relativistic particle accel- 
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erates and pair production occurs. We quantize the fields and use the in/out formalism to 
calculate the pair production, and we introduce the logarithmic negativity [3] to calculate 
the entanglements in different bipartite systems. In Sec. IIVI we consider the case when one 
or both Dirac particles are accelerated. Entanglements between different degrees of freedom 
are calculated and entanglement transfer to the antiparticles will be shown. In Sec. EJ we 
calculate the particle spectrum for scalar particles, and the results are compared with those 
observed by an accelerating observer. We also repeat the calculation of entanglements in 
Sec. |IV] but for scalar particles. A summary and discussion of results and further work is 
given in Sec. EH 

II. ACCELERATING SCHRODINGER PARTICLES 

A free non-relativistic particle with mass m represented by a gaussian wave packet 



^oM) = )= 1 L= exp (1) 



y/4h + 2it/m \ 4b + 2it/m 
satisfies the Schrodinger equation (h = c = 1), 



.cfy = -ld 2 ip 
dt 2m dx 2 

It is reasonable to assume that the center of the wave packet follows the classical traiec- 

n 

tory of the corresponding particle [14j. Therefore, when a linear potential is added to the 
Schrodinger equation, 



.dip -ld 2 i> 

l m = 2^ - max ^ (3) 

we use an ansatz of the form 

V>(x, t) = -x - v t - l -at 2 , t)e iS ^\ (4) 
which when substituted into Eq. (j3J) produces an accelerating wave packet with 



— S(x, t: vq, a) = vqX + axt — -av$t 2 — -a 2 t 3 — -v 2 t, (5) 
m 2 6 2 
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FIG. 1: a) A 3D plot of the absolute value of the two-body wave function fy + Eq. ([6]) at t = 0, for 
v\ = -1, v 2 = 1, x% = X2 = 0, mi = "i2 = Ij cl\ = -0.5, a 2 = 0.5. b) Same as Fig. []Ji, but for £ = 
15. 

where t> o and a are treated as parameters that represent the initial velocity and acceleration 
in the classical limit. 

We write down the two-body entangled wave functions as follows, 



^±{x,y,t) = N[i)(x,t)Vi,a 1 )i)(y,t;v 2 ,a 2 )±i}{x,t]V2,a 1 )'4)(y,t;v 1 ,a2)), (6) 

where N is the normalization factor. For simplicity, we have just set the masses, initial 
positions and widths of the two particles in the wave packet to be the same. If a x or a 2 
7^ 0, the two-body wave function accelerates in a specific direction (see Figs. [T] and [2] for 
examples) . 

We can calculate the purity P of the wave function 



^(x, y, t)*(x', y', t)V*(x, y', t)V*(x', y, t)dxdx' dydy 1 , (7) 



y 
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FIG. 2: Same as Fig. [lb, but for a\ = a<i = (free evolution). 

and hence the Schmidt number: K = 1/P. For a product state, K = 1, and for entangled 
states, K is greater than 1. The entanglement depends on the relative velocity of the two 
particles v instead of the velocity of each of the particles. Therefore, we can choose to 
calculate the entanglement in the frame that v i = and v 2 = v. The Schmidt number for 
ty + (x,y,t) is calculated to be (see Fig. [3]) 



A+ " l + 4//(l + /)" W 

where / = exp(— v 2 ) and v = vmy/b. In this case, the entanglement depends only on the 
product bv 2 m 2 and is independent of the acceleration. This result can be easily understood. 
In the wave function, a appears always as a product with t. As the Hamiltonian is Hermi- 
tian, the evolution operator is unitary, and the entanglement is unchanged under a unitary 
transformation. Thus we can choose to calculate the entanglement at time t = 0, and all 
the acceleration terms will disappear. For K + is equal to 1 when v is equal to 0, and 
K + is equal to 2 when v tends to infinity. For K_ is always equal to 2. Although 
the entanglement is independent of acceleration in both cases, there is a difference between 
the entanglements of ty + and \E'_. The entanglement depends on the orthogonality of the 
two terms in ty±. For a greater v increases the orthogonality of the two terms in the 
wavefunction. The minus sign in cancels the overlap region between the two terms in 
Eq. ([S]), and the remaining parts are orthogonal. As a result, the entanglement of \P_ is 
always maximum in the system. 



K 




FIG. 3: The Schmidt numbers K± as function of the relative velocity of the two particles, made 
dimensionless by v = vm^/b, of the two-body wave functions 

III. RELATIVISTIC FORMALISM 
A. Quantization of Fields 



In order to accelerate a relativistic particle, the Klein-Gordon (KG) or Dirac equation 
with an electric field is considered. A strong electric field makes the vacuum unstable 
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and leads to pair production 
gauge Q] , in Rindler coordinates 
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which has been studied in the time dependent 



201 ] and in a finite region 
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22 



The pair production of scalar particles can be understood in the picture of a wave packet 



12 



241 ] . A wave packet uniformly accelerates from the far past, and then tunneling occurs 
in the region when it meets the potential barrier. The transmission wave packet represents 
the antiparticles while the reflected wave packet represents the particles. Therefore, pair 
production occurs in the tunneling region, and this solves the Klein paradox [25( . Although 
the wave packet formalism is more intuitive, it is not clear how to construct a two-body 
entangled probability density in the Klein-Gordon field. Thus, we quantize the field and 
calculate the Bogoliubov coefficients in the in/out formalism. 



1. Scalar Particles 



The Klein-Gordon equation 
electric field E is 



12 



24| for a unit-charged particle with mass m in a uniform 
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{D^ + m 2 



0. 



(9) 



where = + iA^ and the gauge is chosen to be A = —Ex and A x = 0. We assume 
<f> u (t,x) = Cexp (iut)xuj{x) , and the spatial part solutions are parabolic cylinder functions, 
D_ a _i(x). Details of obtaining the solutions are shown in Appendix [K\ Then we classify 
the solutions in the in/out basis [12j, |26|, such that we have two complete bases to quantize 
the field as 



= E( fl 2 , & + ^ t O ( 10 ) 

U! 

or 

*=£(c*#s ( n ) 

OJ 

where the subscripts p and a label the particles and antiparticles respectively. The operators 
a tjT(&w ) an d o^ip^ ) are the annihilation (creation) operators in the in-basis and out-basis, 
and they are related by the Bogoliubov transformation, 



<C = a*<C-(3*b^, (12) 

where 



27re- i7r / 4 e-^ 2 / 2 







r(i/2 + z/i 2 ) 

e i7r/2 e -7r A t 2 



(13) 



with /i 2 = m 2 /2E, and 



ld 2 - 



(14) 



We can express the in-vacuum state as the linear combination of out states 12] as, 

(3*' 
a 



\0)< 



exp 



a 



|o> 



out- 



(15) 
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We let a = e 1 ^ 1 coshr and (3 = e 1 ^ 2 sinhr, where < r < 0.88 is a parameter related to 
acceleration, and we neglect the phase factors which do not affect the following calculations 
of entanglement. Taking the single- mode approximation, we get the in- vacuum state in 
terms of the out states, 



\0 P )in = CQshr ^tanh"r|n p ) OMf |ra a ) OMt . (16) 

n=0 

Similarly the one-particle state is 

|ljp)<n = To- y^tanh n r\/n + l|(ra + l) P )out\n a )out- (17) 



cosh r 

n=0 

2. Fermions 

For a unit-charged fermion with mass m coupled to an uniform electric field, 

- - rr# = 0, (18) 

where is the vector potential and 7 M is the gamma matrix. Eq. (fl8l) can be reduced to 
two Klein-Gordon equations which are shown in Appendix [Bl The solutions are still the 
parabolic cylinder functions. The in/out basis solution of the second order ODE is still the 
in/out basis solution of the Dirac equation Eq. ffT8l) . Therefore, we obtain the Bogoliubov 



coefficients, which have been calculated in Ref. 



3, 



a n = a f a n ~ Pf° n , I 19 ) 

b n - p f a n + a f b n ', 

where 



,2 



(3f = e-*" , (20) 



01 f = ~ % 




T(i/x 2 

with a j and (3f having the relation, 
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l«/| 2 + l/3/| 2 = l- (21) 
We let atf = cosrye l<?i and /?/ = sinr^, rj being a parameter with values between and 71/ 2 



and related to the acceleration. Also, we can re 
states as in the case of an accelerating detector 



ate the incoming states with the outgoing 



\0 P )in = cosr/e ^|0 p ) cmt |0 a ) out - sin r / |l p ) oui |l a ) ouf , (22) 

|lp)in I ^-p)out 1 0<j)out • 



B. Logarithmic Negativities 



The entanglement can be quantified by the logarithmic negativity 
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29j . For a density 



operator pa,b corresponding to a bipartite system A and B, we define the trace norm | \pa,b\ 
= tr\pA,B\ = tr\ PabPa,b and the negativity 



IIo Ta II - 1 

N e = " P I - , (23) 

where p 7 - 4 is the partial transpose of pa,b with respect to the party A. N e can be calculated 
from the absolute value of the sum of the negative eigenvalues of p Tjv . Then the logarithmic 
negativity of the bipartite system A and B is defined by, 

LN{p AtB ) =log 2 ||2AT e + l||. (24) 
For a product state, LN(pa,b) = 0, and for entangled states, LN(pa,b) > 0. 

IV. ACCELERATING FERMIONS 



Initially, we have the incoming entangled state, 

1 

71 



r— [| 0s,p)m I ^uj,p)in | ls,p)m | lc<j,p)m] • (25) 
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Then either one or both of the particles in u and s modes are accelerated by the electric 
field, and the in states in Eq. ff25]) are replaced by the out states as in Eq. (1221) . If only the 
uj mode is accelerated, we have 



= — ^ {\O s , P )out ® [cosr/e ^10^)^10^)^ - sin r/ll^p)^) l w , a )out] (26) 

"I - 1 ls,p)ottt ® ( I laJ,p) out | ^ui,a) out ) J • 

If both the s and uj modes are accelerated with the same ry, we have 



= -/={[cosr/e 4<?il |0 Si p) out |0 Sia ) o ^ - sinr / |l SiP ) OU i|l Sia ) out ] (27) 
(8) [cos r/e"^ 1 1 WiP ) outfit ,a) out - sinr/ll^o^l^a)^] 

+ \{\^s,p) out\^s,a) out) ® ( | ^-w,p) out \ ®u,a) out))} ■ 

The degradation of entanglement in the case of an accelerating detector is due to the fact 
that some degrees of freedom have been traced out. An accelerating detector 'sees' the space- 
time being split into two causally disconnected regions, and it cannot access information in 
one of them. We have verified explicitly that the entanglement between the particles in s 
mode and uj mode is unchanged if there is no tracing out of any space-time region. On the 
other hand, in the case of accelerating particles, the detector, which is in an inertial frame, 
can access all degrees of freedom and the orthogonality of the states is unchanged; therefore, 
the entanglement of accelerating particles is unchanged. 

However, more degrees of freedom are produced and we can calculate the entanglements 
between different bipartite systems. In Ref. [4|], it was shown that entanglement is Lorentz 
invariant. If one traces out the momentum, the entanglement decreases, and the entan- 
glement is transferred from the momentum to the spin degrees of freedom. We will show 
that entanglement transfer also occurs in accelerating fermions, from the particles to the 
produced antiparticles. 

If only the particle in the uj mode is accelerated, we can study the three bipartite systems: 
A = the s mode, B = the particles in the w mode, the antiparticles in w mode, or the entire 
w mode including both the particles and antiparticles. The density matrices are called p SjP , 
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p Sja , and p s ,(p,a) respectively. The entanglements are 



LN(p s> (p, a )) = 1, 



< LN(p s , p ) =log 2 (l + cos 2 r / ) 
LN{p s>a ) =log 2 (l + siii 2 r / ) 



(28) 



v 



which are plotted in Fig. HI It is obvious that the entanglement of p SiP is transferred to p Sja . 

When both the s and us modes are accelerated with the same 77, we can calculate the 
entanglements between the five bipartite systems: particles in s mode and particles in u 
mode (pp, p ), antiparticles in s and antiparticles in u (p a , a ), antiparticles in s and particles in 
uj (p a ,p), particles in s and antiparticles in u (p p , ), an d the entire s and u> modes (/O(p,o),(p,a))- 
The logarithmic negativities are 



By symmetry, LN(p a ^ p ) = LN(p p ^ a ). The results are shown in Fig. HI The entanglement 
is transferred from p pp not only to p p>a , but also to p a<a . In fact, when the acceleration 
of the particles tends to infinity, the entanglement is completely transferred to between 
the antiparticles p a ^ a . Note that in both cases, the negativities of the subsystems add 
up to the that of the total system, i.e., N e [p SiM ] = iV" e [p SiP ] + N e \p 8>a ], Ne[p( Pj a),(p,a)} = 



A. Spectrum 

The relation between the in states and out states for scalar particles is just the same 

n q 

as the relation between the Minkowski states and Rindler states in p, LDj . However, their 
spectra are different. For both cases of accelerating particles with an inertial observer and 
inertial particles with an accelerating observer, the spectra are 



LN(p (yPya) ^ a) ) = 1, 

LN(p PtP ) = log 2 [1 + cos 4 77] , 
LN(p ata ) = log 2 [l + sin 4 77] , 
k LN(p Pia ) = log 2 [l + cos 2 Tf sin 2 77] . 



(29) 



N e [p p J + N e [p a<a ] + N e [p p<a ] + N e [p a>p }. 



V. ACCELERATING SCALAR PARTICLES 



(oicxrio) 




(30) 
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FIG. 4: Logarithmic negativities of several bipartite systems when one or both fermions are accel- 
erated, the magnitude of which is parameterized by ry. In both cases, the entanglement between 
the entire s mode and uj mode is unchanged (dot-dashed line). The solid lines show the results 
when both particles are accelerated together, for three bipartite systems: particles in s mode and 
particles in uj mode (p P , P ), particles in s mode and antiparticles in uj mode (p P , a ), and antiparticles 
in s mode and antiparticles in uj mode (p a ,a)- For comparison, the dashed lines show the results 
when only the particle in the uj mode is accelerated, in which case the two bipartite systems are 
particle in s mode and particles in uj (p s , P ), and particle in s mode and antiparticles in uj {p s ,a)- 

The spectrum of accelerating particles with an inertial observer is S w = exp(— rem/ a), where 
a = E/m corresponds to the acceleration of the particle in the classical limit. However, 
a uniformly accelerating detector measures a spectrum S u = l/[exp(27r|c<j|/a) — 1]. In the 
classical limit, the detectors observe the same particle trajectories; however, a uniformly 
accelerating detector measures a different spectrum of particles as that by an inertial detector 
on uniformly accelerating particles. 

B. Entanglement 

For scalar particles, we have the same initially entangled state in Eq. fl25|) . If only the 
particle in uj mode is put in a uniform electric field, the entangled state becomes 
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r 



0.2 0.4 0.6 0.8 

FIG. 5: Same as Fig. [H but for scalar particles. Again, the entanglement of the entire s mode and 
uj mode, indicated by the dot-dashed line, remains unchanged. Note that LN(p s>a ) = LN{p pA ) = 
LN(p^ a ) = for all r. 



*/ = ^ { K,s)out 



1 00 

~r — ^ ] tanh r\n PiUJ ) out\ n a,u>) out 



coshr 



n=0 



+ 



(31) 



I lp,s) 



out 



I 00 1 ] 

-^y^t&n\i n ryJn + l\(n + l) p ^) out \n a ^) out \ . 

cosh r 

n=0 J ) 



If both particles in uj and s modes are put in the electric field with same r, we have, 
1 



tanh™ r |n PiS ) ou4 \n a , 3 ) 



out 



n=0 



cosh 2 r 

oo 

tanh™ rVra + 1 1 (n + l) p , s ) out\n a , s ) 



tanh™ r\n P:LU ) out \n a ^) 



out 



n=0 



(32) 



n=0 



cosh 4 r 

oo 

^ t&nh n rVnTl\(n + l) p>ul )out\n a ,u>) 



out 



out 



n=0 



More degrees of freedom have arisen due to pair production. Again, the entanglement 
between the entire s and lj modes remains unchanged, such that LN(p S)U1 ) = LN{p s ^ j0 ^) = 
LN(p(p ia ) i (p ja )) = 1 for all r. This is because the Bogoliubov transformation is linear, and 
the orthogonality property of the states is unchanged. 

We calculate the entanglements of p S}P , p s<a , p PtP , p a<a and p Pja . To calculate the density 
matrix, p SjP (p s , a ), we trace over the antiparticles (particles). We take the partial transpose, 
Ps P i by interchanging the s mode's qubits to get an infinite block-diagonal matrix. The (n, 
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n + 1) block matrix is, 

Tann r I sinh 2 r coshr 

2 cosh 2 r 1 v^+T tanh 2 r 

y cosh r 

Then we calculate the negative eigenvalues from each block matrix and obtain the logarith- 
mic negativity of p SiP , 




LN (Ps, P ) = log 2 

We also calculate the LN(p P}P ) numerically. The results are shown in Fig. [5j 

In contrast to fermions, there is no entanglement transfer to the antiparticles for scalar 
particles, and LN(p Sja ) = LN(p Pia ) = LN(p a}a ) = for all r, even though the entanglement 
between the particles in the s and to modes decreases as r increases. 



1 

2 cosh 2 r 



E 

n=0 



tanh 2n r 



2 cosh r V \ sinh r 



+ tanh r + 



4 

cosh 2 r 



(34) 



C. Entanglements if the Number of Produced Pairs is Restricted 

However, if we constrain the number of pairs produced, LN(p s a ), LN(p p a ) and LN(p a ^ a ) 
are all nonzero. If only M pairs can be produced in a mode, 



\0)in = 7— y^t&nh n r\n p ) out \n a ) out , (35) 

coshr ^— ' 

n=0 

|l P )m = — t%- ^ tanh" ry/n + l|(n + l) p )out\n a )ouu 
cosh r ^— 

n=0 



where Ni and N2 are normalization factors, 



N 1 = (l-tanh 2M+2 r)" 1/2 , (36) 
N 2 = [i _ ( M + 1) tanh 2M r + M tanh 2M+2 r] ~ 1/2 . 

We take the partial transpose of p SjP which has M diagonal block matrices and the nth 
block is, 



tanh 2 - 4 r / ^(n - 1) ™f^v^ 

2 cosh 2 r I ^^tanhV Ar 2 tanh 4 r 

\ coshr v i 
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(37) 



LN 




0.2 0.4 0.6 0.8 



FIG. 6: Same as Fig. \5\ but for M = 1 (solid lines) and M = 2 (dashed lines), and p s>a is shown 
as well. 

We then sum up the negative eigenvalues of the nth blocks to calculate the logarithmic 
negativity 



LN{ Ps 



log 2 



M 

£ 

n=l 



tanh 



2n-4 



2 cosh r 



(n - l)iV 2 2 
cosh 2 r 



+ JVi tanh 4 r- 



(38) 



(n ~! 2 )jV| +^ 2 tanh 4 r 
cosh r 



+ 



±NlN% tanh 4 r 
cosh 2 r 



The partial transpose of p s>a also has a block diagonal structure and only the last block 

tanh 



N1N2 tanhr 



cosh ) 



2 COsh r \ AtAjimii. >• .77 Q 

contributes to the negative eigenvalue. Then the logarithmic negativity of p s>a is 

iVftanh 21 ' 1 - 2 ' 



(39) 



LAT(p Sia ) = log 2 <^ 1 - 



1- Wi + 



4iV 2 M tanh 2 r 



A^ 2 cosh r 



(40) 



2 cosh r 

We show the results M = 1 and 2 in Fig. [6j In both cases, LN(p s ^ p ^) is equal to 1. 

When M is finite, the entanglement of p s>a is not zero and increases with r while that of 
p S)P decreases with r. However, the entanglements of both p SjP and p SA are reduced if more 
particles are produced (M increases). The dependence of LN(p s ^ p ) and LiV(p s a ) on M at 
infinite acceleration are shown in Fig. [71 Although the entanglement of p SyP decreases more 
for greater M, the entanglement of p Sya also decreases with a similar trend and goes to zero 
when M — > 00. Therefore, there is no transfer of entanglement to the antiparticles for scalar 
particles when the number of produced pairs is not restricted. 
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FIG. 7: Logarithmic negativity of p StP and /9<j ja at infinite acceleration (r = 0.88) as a function of 
the number of produced pairs M. 

We next calculate the case when both particles in uj and s modes are put in the uniform 
electric field such that they are both accelerated and have the same r. In the out basis, the 
state in Eq. ( 1251) becomes, 

M 1 r m 

tanh r\n p! s) ou t\n a ^ s ) ou t 



1 



n=0 



cosh r 



\pl [ cosh 

M 

y^tan^rVra + l\(n + l) PtS ) ut\n a , s ) 



out 



n=0 



+ 



22 tanh™ r | n v ^ ) out \ n a ^ ) f 

n=0 
M 

J^tanh n r\/n + l|(n + l) P ,o;) OU i|^a^) 



n=0 



Now, we can consider even more bipartite systems between the two modes. We calculate 
the entanglements between the particles or antiparticles in the s mode and the particles 
or antiparticles in the u mode, i.e., of p PtP , p a>a , P P ,a and p a , p for the cases M — 1 and 
2. As expected from symmetry, LN(p p ^ a ) is equal to LN(p a ^ p ). The results are shown in 
Figs. [8] and [9] for M = 1 and 2 respectively. As in the case when only one of the particles is 
accelerated, the entanglement between the particles in s and u modes is degraded when M 
increases. It seems that entanglement transfer also occurs for scalar particles for finite M. 
However, in this case we note that the negativities of the p a)a , p PA , p a ^ p , and p PjP do not sum 
to a constant, and we cannot identify a causal relation between the decrease in entanglement 
between the particles and the increase in entanglement between the antiparticles. A more 
rigorous definition of entanglement transfer is needed before we can discuss the issue for 
scalar particles further. 
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0.2 0.4 0.6 0.8 

FIG. 8: Logarithmic negativities of different bipartite systems when the particles in both u> and 
s modes are accelerated, indicated by p PtP (between the particles and particles), p PA (between 
particles in s and antiparticles in u> modes), and p aj(l (between antiparticles in s and uj modes), 
for M = 1. For comparison, the previous results (shown in Fig. [6|) when only the particle in u is 
accelerated are plotted as dashed lines. 




FIG. 9: Same as Fig. El but for M = 2. 

VI. CONCLUSION 

We have studied how the entanglement of a pair of particles is affected when one or both 
of the pair is uniformly accelerated, as measured by an inertial detector, and compared the 
results with that of inertial particles observed by a uniformly accelerating detector. While 
there is a degradation of entanglement in the latter case due to the splitting of the space- 
time, the entanglement in the former case is unchanged by the acceleration when all degrees 
of freedom are considered. Furthermore, the spectrum of the uniformly accelerating particles 
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is different from that seen by uniformly accelerating detectors. 

When one of the particles - the one in u> mode - is uniformly accelerated, the entanglement 
is transferred to the produced antiparticles for fermions, while there is no such entanglement 
transfer for scalar particles. For scalar particles, when the number of produced pairs is 
restricted, the entanglement of p s>a increases with r. However, at any r, LN(p sa ) decreases 
as M increases and goes to zero when M — > oo. 

When both particles in s and u modes are uniformly accelerated (with the same r or 
rf), we have even more bipartite systems. For fermions, p a ^ a takes up all the entanglement 
at large acceleration. For scalar particles with restricted number of produced pairs, the 
entanglements of p PA and p a ^ a increase with r. However, if there is no restriction of the 
number of produced pairs, no entanglement transfer to the antiparticles is observed for 
scalar particles. 

Our results raise the possibility that when an entangled pair falls into a black hole, their 
entanglement may be partially transferred to the produced particles, which should not be 
ignored in considering the black hole information paradox. Studying quantum entangle- 
ment in curved space-time may therefore give us insights on the relation between quantum 
mechanics and general relativity. 

APPENDIX A: THE BOGOLIUBOV COEFFICIENTS IN THE SCALAR CASE 

We assume the form of solution of Eq. flU]) as 



(Al) 



where C is a normalization constant, and we obtain from Eq. 



dx 2 



+ E 2 (x - tu/Ef 



The solutions of Eq. (]A2I) can be found in |30j , and they are parabolic cylinder functions 



(A2) 



D^[±(l+i)y/E(x-u)/E\, 
D_^_Ml-i)VE{x-uj)/E]. 



(A3) 
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We can use the saddle point method to classify the solutions in the in/out basis [l2j, |26] 
and have the in-basis functions, 

-37i7i 2 /4 

(2£)W 



t) = e -—^e^D^_ 1/2 [e-^V2E(x - u/E)], (A4) 



<j>l n a (x,t)=^ p (-x 1 t), (A5) 

where /i 2 = m 2 /2E. The subscripts p and a stand for particles and antiparticles respectively. 
We also obtain the out-basis solutions, 



CM) = r; P (-*,-*)- ( A? ) 

The solutions have been normalized by the Klein-Gordon scalar product, 

J dx(f>%* p (x, t)W t (f)™ p (x, t) = +5(co - J), (A8) 

J ^;(m)Sc(m) = -s{u-j). 

As there are two different complete bases, we can quantize the field in two ways, 



= Ek>^+ 6 ^c) 



(A9) 



or 



OJ 

From [30[, we can get a relation between the out-basis and in-basis solutions, 



(A10) 



D + i)VE(x-u/E)} 



(All) 



£-i-^[(l + 0v^(x-aV£)] + 



27re^ i7r / 4 e-^ 2 / 2 



^_i + iu2[— (1 — i)\fE{x — L)/ E)\. 



We write Eq. (lAllj) in the form of in-basis and out-basis solutions explicitly, 



b out = 3(b in * + acb in . 



(A12) 



Since (3 = e 1 ^ 2 sinhr, r = sinh 1 [exp(— nm/2a) 
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APPENDIX B: REDUCING THE DIRAC EQUATIONS TO TWO KLEIN- 
GORDON EQUATIONS 



From Eq. (1181) . we let 



(Bl) 



to obtain 



[(p _ A) 2 -m 2 - ia^( Plx - A„)(p v - A v )]<j> = 0, 



(B2) 



where o^ v = zpy^, 7"]/2. We choose the gauge to be A = 0, A 3 = —Et. We then substitute 
the potential in Eq. flB2j) to get 

)■> / c \ 2 



0. 



where 



«3 



a 3 
a 3 



and (T3 is the Pauli matrix. We then assume that the solution has the form 



(B3) 



(B4) 



(f) = e lkz n(t) 



(B5) 



where 



with the spinors, 



Ul 



u 2 



«3 



«(*)=X}/x(t)UA, 
A=l 



J_ 
1 

w 
( 1 ^ 

1 o 
-i 

V o / 



V2 



,U4 



1 

71 



/ o \ 

1 



V- 1 / 
/o\ 
i 



V 1 / 



(B6) 



(B7) 
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With the relations, 



a 3 u x = 7]u x , 



we can get two Klein-Gordon equations, 



1, for A = 1,2, 
-1, for A = 3,4, 



(B8) 



dt< 



+ E 2 [t 



k 

E 



+ m 2 — %Er\ 



fx(t) = 0. 



(B9) 



The solutions are parabolic cylinder functions. 



We can then classify the in/out solutions [24j, |3l|. / 3 and / 4 are dependent on fi and / 2 , 



and so we just consider the cases of A = 1, 2 in the following. Neglecting the normaliazation 
factors, we write down the in/out solutions, 



km 
V 



e ikz D llA 2[-(l-tWE(t-k/E)}, 
in = e «* jD _. Ma _ 1 [_(i + i)y/E(t - k/E)}, 
out = e i**£,_. Ma _ 1 [(i + _ k/E)], 

e ik *D ilx2 [(l-i)^/E(t-k/E)\, 



(BIO) 



where fi 2 



nr 



l /2E. We substitute the solutions in Eq. (1B10I) into Eq. ( 1B1I) . to obtain the 



solutions in the Dirac equation, i.e., Eq. ( flBl) . We show the calculation of $™ with A = 1 as 
follows. 



% n i = [Yip^ - A^) + m} Ul e lkz D^[~(l - t)VE(t - k/E)} 



J.kz 



(Bll) 

muiD ifl 2[-(l -i)VE(t - k/E)] - (1 - ^/iV-E^A^-ihil -i)VE(t - k/E)} . 



where u l = j Ui = 7 U\. We normalize it and calculate "0* n , and ip° . We can relate 
the solutions using two mathematical relations 



[l+i)x], (B12) 



2 \/2ne~ nlJ - 2 / 2 



r(i-z^ 



A^[(l - *>] = e"^ A^h(l - *>] + ^ . 2/n , £>-i-^[-(l + i)x], 



r(-^ 2 /2) 



and so we can obtain the Bogoliubov coefficients, 



■TTfJr 



(B13) 



/x 2 T(ifj, 2 ) 
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Since f3f = sinr/, r/ = arcsin [exp(— irm/2a)]. 
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